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^ ! Abstract 

H 

Cd ' We postulate that the Fermi function should be derived from the amplitude, not from the solution of 

the Dirac equation, in the quantum field theory. Then, we obtain the following results. 1, We give the 

amplitude and the width of the neutron beta decay, n ^ p + e~ + v^ to one loop order. It is carried out 

by the Feynman parameter integration. 2, As the result, we find the terms which can be interpreted as the 

Fermi function expanded to order a. 3, We also give the same result using complex analysis. 4, We check 

that there are no such terms in the similar process, i/g + P — >■ e"*" + n. 5, We perform the Fermi function 

"■^ ■ expanded to order a^ using complex analysis. 
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The calculation of beta decay rates, for example n -^ p + e + i/g , requires the Fermi function. The Fermi 
function represents the effect that the electron runs through the electromagnetic potential caused by the 
proton. This function affects the beta spectrum, the decay width, and the lifetime of the parent particle. 
The Fermi function has been derived so far as the solution of the Dirac equation in electromagnetic 
VO ' potential [l| [4I |3| [J] . In beta decay, the decay itself is caused by the weak interaction and treated as the 

\l . intermediate state, which is represented as the amplitude^ The final state particles are in electromagnetic 

^^ ' potential. This effect is factorized as the Fermi function [5]. For the non-relativistic limit in neutron beta 

CO . decay, it takes the form 

o 

(N : _ 2TTa/v 



where a is the fine structure constant and v is the electron velocity in the neutron rest frame. Calculating 
^ . the decay width, the Fermi function is multiplied by the absolute square of the amplitude, and integrated 

over the phase space. This is the same for the loop amplitude. For the sake of simplicity, in this paper, we 
set the parent and daughter nucleons as neutron and proton, respectively. 

From a quantum field theoretical point of view, this potential effect is also the interaction of intermediate 
state. Generally, the created particles should be considered as the asymptotic fields in the far future. 
The asymptotic fields are free from the interactions between the created particles. The electromagnetic 
interaction should be derived from the loop diagrams. Furthermore, this effect should appear in perturbation 
theory as 
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order by order with respect to a. We note that F^ji is actually expanded with respect to a/v rather than 
a. 

The tree level beta decay diagram is 

p+{p') 



e-{q') =iMo, 

i>e{q) 
where the parameters in each parenthesis represent their momenta. Also, the one loop diagrams are 

p+{p') 

n{p) " --^ >-^ e-(q') ^ iMiL 





(3) 

and the field strength renormalization. According to Ref. [6|, the one loop amplitude is divided into three 
parts as iAfiL — iMi + iM2 + iM^. We calculate only iAIi + iM2, which do not depends on a^'^ki, in the 
numerator of proton propagator, where cr^'^ = i^j^j''' — '^'^j^)/2. To cancel the infrared divergence, we 
consider the sum of two bremsstrahlung diagrams, iMh. 

The detailed calculations are given in Appendix [Appendix A[ and the result is 
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where Gp, w„, m^, me, and E^ are the Fermi constant, neutron mass, proton mass, electron mass, and 
electron energy, respectively; fcjv/ = to„ — mp — E^ represents the maximum photon energy for given E^ ; C 
represents the Gamow- Teller coupling constant relative to the Fermi transition; Li{x) is the Spence function 
defined as 



Li{x) 



dz 



log(l - z) 
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We set nin — mp ^ Ee,me and the neutrino is massless. dV^ is a part depending on iM^, which does not 
depend on i) [6|. We do not calculate dT^ in this paper. The electron velocity is represented as z; = \q'\/q'o 
in the neutron rest frame. The first term in the curly brackets of Eq. Q is inversely proportional to v. This 
term is interpreted as the one corresponds to order a Fermi function expressed in Eq. ^. We note here, 
the Eq. Q differs 1/v term and the constant from Ref. [6|. 

According to Appendix [Appendix A[ the amplitude is approximately written as 
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for w <C 1, where /i is the photon mass introduced to regulate the infrared divergence. The last term does 
not affect the one loop width, but affects the two loop one as explained later. 

Here, we extract Eq. (|6]) using complex analysis. According to Appendix [Appendix A /5a contains 
this term. I^a originates from one loop diagram in Eq. (jH]), not from the bremsstrahlung or field strength 



renormalization terms. Furthermore, I^a does not contain k in its numerator. We explain this reason later. 
Therefore, we start from 

r d^k -Ae^Mpp'-q' 1 1 

* ^^y (27r)4 (p' _ fc)2 _ ^2 + ie (q' + k)^ - ml + ie k"^ - ^i"^ + ie' ^' 

where ie represents the Feynman prescription; e is the electromagnetic coupling constant. 

Next, we integrate Eq. ([7]) with respect to /cq, where the contour integral is taken clockwise. The locations 

of poles are fco = Pq + Jp'^ +k? — ie, fco = — Qq + V 9o + k"^ — 2<7' ■ k — ie, and fco = \/k^+Ji^—ie. Since 
we focus on the terms proportional to a/v, which dominate for v ~ \q'\/q'a <^ I, we set p'q ^ q'q ^ \q'\. The 
contributions which converge for i; ^ can be ignored. For |fc| > q'p, the integral over |fc| converges even if 
we set u — >■ 0. Hence, we set ^q ^ I '''I: ^nd the locations of poles are fco — 2pg, fco ~ j-r{k^ — 2q' ■ k — ie), 

and fco — \/k'^ + /i^ — ie. Writing only the contribution from the second pole explicitly, 

/cfik 1 1 
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where the part, "+••." represents the contribution from other poles. 

Defining x — k/\q'\, \x\ = x, cos 9 = q' ■ k / {\q' \\k\) , and p, — n/\q'\, the integrand becomes dimensionless 
as 

, ^ e'^iMo q'n f dcosddx x^ aiMo f°° , logfa;^ + 2x ~ ie) - logCx^ - 2x - ie) 
iMiB——r^^— : — - + ---3 xdx-^ ^ ^^ ^ 
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In the right hand side, the part " + •••" is ignored since these terms are not proportional to the factor 1/v. 
Generally, the dimensionless integrand does not contain the factor v. Then, it can be ignored that the terms 
which do not have the factor 1/v as a coefficient of the integral. Here, if the numerator of the integrand 
have the term which contains fc, it gives additional factor \q'\ when we take the integrand dimensionless. 
Then, such a term cannot be the candidate of the Fermi function. 

We use complex analysis to the x integral. Since the integrand in Eq. ^ is the even function, the 
amplitude can be written as 

iMi 3-—iMo / ^5 r:^ [log(a;2 + 2x ~ ie) - log(a;2 - 2x ~ ie)] . (10) 

27rw J^oo x"^ + fi'' ~ le ^ ■' 

To apply the residue theorem, we integrate Eq. (I10[) by parts to form 

a f°° x"^ 

iMi 3—iMo / dx^^—- . log(x2 +p.^- ze). (11) 
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Applying the residue theorem, the amplitude becomes 
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where the contour is depicted in Fig. [TJ 

This result is consistent with Eq. ^ not only the real part but also the imaginary part in the curly 
brackets. 

More than one charged particles does not exist at the same time in the scattering process Pg +P -^ e^ + n. 
Therefore, the electron and proton are not affected by the electromagnetic potential and the amplitude should 
not contain the a/v terms which can be interpreted as a part of the Fermi function. Here, we verify it. 

According to Ref. [7|, the diagram is 




Figure 1: x contour and the location of the poles 
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Extracting the related terms by the similar manner as the beta decay, 
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We integrate with respect to k^ along the counter-clockwise. After the same approximation with the beta 
decay, the poles we focus on are k^ ~ ~2p^(^^ ~ *£) ^md fco = ~'^{^'^ ~ 2?' -k — it). Each residue is the 
same value with the opposite sign. Then, the amplitude becomes 
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where iM^ is the tree level amplitude. 

As a result, this scattering process does not have a/v term, which can be interpreted as the part of the 
Fermi function. This is because the signs on k in the electron and proton propagators are the same. It is 
just equivalent to the two charged particles do not exist at the same time. This supports that the a/v terms 
in beta decay correspond to the potential effect. 

Two loop ladder diagram is 
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where fci and ^2 are the photon momenta, respectively. There are some other two loop diagrams. However, 
we are only interested in the term proportional to (a/v)'^. This term originates only from iM2L- After the 
similar calculation to the one loop amplitude, the amplitude contains 
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Summing this equation, Eq. (J12p . and the tree level amplitude, we give 

zMo + zAfiL + iM^L 9*Afo 1 + log ^ - — -— TT^ + - log2 ^ U . (15) 
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The absolute square of them is 

|zMo + zMii + iMsLi' 9 \iMof (l + — + ^) + Oia'). (16) 

The logarithmic terms in Eq. (fT5)) are canceled. For w <C 1, the decay width has the form 

dr-dr3(xi + — + ^r^- (17) 

This is consistent with the Fermi function up to order a^ 
We conclude the main results as follows. 

1. We reviewed the one loop beta decay amplitude to find the terms proportional to a/ v. It can be 
interpreted as the part of the Fermi function. 

2. The scattering process Ve + P ^ e^ + n does not have such terms. 

3. We give the result that the two loop beta decay amplitude has the terms proportional to (a/u)^. These 
are consistent with the expanded Fermi function up to order a^ . 

These results suggest that the potential effect named the Fermi function should be considered as a part 
of the amplitude. 

To confirm our conclusion, it is necessary to carry out the Fermi function to higher order. If the systematic 
calculation will be carried out, we can sum up all the order of contributions. 

For a/v > 1, the perturbation up to the finite order does not work. We must sum up all order of a/ v. 
The result should become the full Fermi function written in Eq. ([T]). Then, we propose the decay width to 
form 
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In two loop calculation, we were only interested in the terms proportional to (a/v)"^. However, the terms 
proportional to a^/v may exist. These terms also affect the decay width. We should consider them for the 
higher order calculation. 

Our result does not affect the practical use except for v < a. For instance, this result only slightly affects 
the Kurie plot [8| and Ref. [9|. However, the theoretical calculation of the nuclear lifetime is changed. 

This study can be applied to beta decay of other nuclear species by exchanging a — ?> Za, where Z is 
the atomic number of the parent particle. The loop diagrams which contain the photon propagator between 
parent nuclear and the daughter particles do not give the a/v term as explained in the calculation of the 
scattering process. Also, we confirmed that the {a/v)^ term does not appear in the corresponding two loop 
diagrams. Our study is more important for the beta decay of nuclear which have larger Z. 

Appendix A. the Detail of One Loop Calculation 

The tree level amplitude is 

iMo = -'-^u{p'){l - Cj')u{p)u{q'){l - j')v{q). (A.l) 



According to Ref. [6], the one loop amplitude can be separate in three parts as 

iMiL = iMi + iM2 + iMs- 
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iMi picks up the factors {2q' + k)^ from the electron propagator and (2p' — k)^ from the proton propagator. 
iM2 picks up the factors a^'^ki, from the electron propagator and {2p' — k)^ from the proton propagator. 
iMa picks up the remaining factors. 
The tree level width is 
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where E^, and to„ are the neutrino energy and the neutron mass, respectively; ^^ represents the spin sums. 
The one loop width is 
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where SZ^ and 5Zp are the one loop electron and proton field strength renormalization, respectively; dVb is 
the term originated from the bremsstrahlung; dT^ is the term originated from iA/3. 
Here, iM2 is written as 
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Then, the cross term between iM^ and iA/2 is 
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Then, the decay width takes the form 
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where li are defined in Appendix [Appendix A.^ 

Appendix A.l. dVb 

The bremsstrahlung amplitude is approximately written as [lO| 
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for small fc, where e^(fc) in the numerator represents the polarization vector of the external photon. 
The absolute square of this amplitude is 
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where j3 = |fc|/A:o and w — k ■ q' /{\k\\q'\). Here, we define k = |fc|, and 
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where jjl — \/k'^ ~ I'^'P is the photon mass. Then, the bremsstrahlung part is 
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Therefore, the Eq. (jA.6p becomes 
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According to Appendix [Appendix A.2[ the one loop width finally takes the form 
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Appendix A. 2. Ii ^ Iq 

Using the Feynman parameter integral, we define /^'s. Here, /i + /s is derived from 5Zp/2. Similarly, 
I2 + Ia is derived from 5Ze/2. Also, /5 + Iq corresponds to iMn,. These are already canceled the ultraviolet 
divergence. The results are as follows. 
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We note here that b/c — 1/v. 
The sum of these li 's is 
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The last two terms diverge for u — > 0. The latter one, which contains /i, does not affect the one loop decay 
width. However, this term has an nontrivial, important roll to two loop calculation. 

References 

[1] E. Fermi, Z. Phys. 88, 161-177 (1934). 

[2] Fred L. Wilson, American Journal of Physics 36, 1150 (1968). 

[3] H. R. Hulme, Proc. R. Soc. London A133, 381-406 (1931). 

[4] Frank. P. Calaprice, and Barry. R. Holstein, Nucl. Phys. A273, 301 (1976). 

[5] D. H. Wilkinson, Nucl. Phys. A337, 474-504 (1982) 

[6] Masataka Fukugita, Takahiro Kubota, Phys. Lett. B598, 67-75 (2004). e-Print: 'hep-ph/0405298| 

[7] Masataka Fukugita, Takahiro Kubota, Acta Phys. Polon. B35, 1687-1732 (2004). e-Print: hep-ph/0403149 

[8] Susan Gardner, Veronique Bernard, Ulf G. Meissner, Phys. Lett. B598, 188-196 (2004). e-Print: hep-ph/0407077] 

[9] Akihiro Matsuzaki, Hidekazu Tanaka, Prog. Theor. Phys. 123, 1003-1012 (2010). e-Print: arXiv:0908.4163r^ 

[10] Toichiro Kinoshita, Alberto Sirlin, Phys. Rev. 113, 1652-1660 (1959). 



